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ABSTRACT 

Since giant planets scatter planetesimals within a few tidal radii of their orbits, the loca- 
tions of existing planetesimal belts indicate regions where giant planet formation failed in 
bygone protostellar disks. Infrared observations of circumstellar dust produced by colliding 
planetesimals are therefore powerful probes of the formation histories of known planets. 
Here we present new Spitzer IRS spectrophotometry of 111 Solar-type stars, including 105 
planet hosts. Our observations reveal 11 debris disks, including two previously undetected 
debris disks orbiting HD 108874 and HD 130322. Combining our 32 /im spectrophotome- 
try with previously published MIPS photometry, we find that the majority of debris disks 
around solar-type stars have temperatures in the range 60 < Td us t ^5 100 K. 

Assuming a dust temperature Td us t — 70 K, which is representative of the nine debris 
disks detected by both IRS and MIPS, debris rings surrounding Sunlike stars orbit between 
15 and 240 AU depending on the mean particle size. Our observations imply that the planets 
detected by radial- velocity searches formed within 240 AU of their parent stars. If any of 
the debris disks studied here have mostly large, blackbody emitting grains, their companion 
giant planets must have formed in a narrow region between the ice line and 15 AU. 



1. INTRODUCTION 

Sometimes the giant planets that failed to form re- 
veal as much about the birth environments of other 
Solar Systems as the giant planets we actually find. 
Our objective is to characterize circumstellar dust 
produced by colliding planetesimals orbiting planet 
hosts and use the orbital radii of such planetesimals 
to constrain the locations where the known giant 
planets may have formed. Radial- velocity searches 
do not yet have the time baseline to detect giant 
planets more than ~ 7 AU from their host stars. 
Nevertheless, we expect that giant planets do not 
reside within observed debris disks — which typically 
orbit at 10 AU or more — because giant planets scat- 
ter planetesimals within a few Roche radii of their 
orbits (e.g. Bryden et al. [2000], Thommes et al. 
[2003], Ida & Lin [2004], Goldreich et al. [2004]). As 
long as planetesimal belts don't migrate on a large 
scale, grain temperature measurements that con- 
strain the location of circumstellar dust highlight 
the regions in long-gone protostellar disks where gi- 
ant planet formation didn't succeed. 

The "Vega Phenomenon" — an excess of infrared 
emission caused by circumstellar dust — was first ob- 
served by Aumann et al. (1984) using the Infrared 
Astronomical Satellite. Since then, observations 
by the Spitzer Space Telescope have revealed the 
temperatures and locations of scores of extrasolar 
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Kuiper Belts (e.g. Chen et al. [2006], Hillenbrand 
et al. [2008]). Here we combine new Spitzer IRS 
spectrophotometry of planet hosts with previously 
published MIPS data, allowing us to constrain grain 
temperatures and planetesimal belt locations and 
exclude regions inhabited by circumstellar debris 
from the possible giant planet formation zone. 

Massive planets with small orbital radii almost 
certainly formed at larger radii, where they had ac- 
cess to abundant solid material (Bodenheimer et 
al. 2000), and migrated inwards through either dy- 
namical interactions with massive gas disks (Lin et 
al. 1996) or Kozai cycles triggered by a companion 
(Libert & Tsiganis 2009). Migration will disrupt 
any planetesimal belt in the giant planet's path, in- 
crease the collision rate of planetesimals, and ac- 
celerate the evolution of the planetesimal belt to 
reduce the debris production at older ages. Indeed, 
migration of the gas-giant planets in our own Solar 
System and the subsequent disruption of the Kuiper 
Belt is one possible explanation for the Late Heavy 
Bombardment that occurred in the inner solar sys- 
tem 3.9 Gyr ago (Gomes et al. 2005). Raymond et 
al. (2006) found that if migration occurs quickly and 
the planet settles at an orbital radius < 0.25 AU, the 
planetesimal disk can regenerate to form terrestrial 
planets, and consequently may form a long-lived de- 
bris disk. However, disk regeneration is suppressed 
if the migration stops at a radius a > 0.5 AU and 
the debris may dissipate rapidly. 

The evidence for brighter debris disks around 
planet hosts than non-hosts has at times been in- 
triguing (e.g. Beichman et al. 2005, Kospal et al. 
2009), but most studies show no statistically sig- 
nificant differences between the circumstellar dust 
populations of planet hosts and other stars (Moro- 
Martfn et al. [2007], Bryden et al. [2009]). Despite 
a few compelling exceptions — for example, the ex- 
istence and location of Fomalhaut b was predicted 
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through analysis of the shape and surface brightness 
profile of Fomalhaut's debris disk (Quillen [2006], 
Kalas et al. [2008]) — statistical analysis indicates 
that known giant planets tend to be dynamically de- 
coupled from any planetesimal belts orbiting their 
host stars. 

Since 60% of the exoplanets discovered at the time 
of this writing have semimajor axes > 0.5 AU 6 , we 
expect that few planetesimal belts in the path of 
migrating planets were able to regenerate accord- 
ing to the Raymond et al. (2006) scenario. The 
simplest hypothesis, then, is that migratory plan- 
ets formed nearer their host stars than the dust- 
producing planetesimal belts generally studied with 
IRAS and Spitzer. Migrating inward, they left the 
planetesimal belts intact. Since giant planet for- 
mation efficiency decreases with distance from the 
star (Safronov [1969], Pollack et al. [1996]), we ex- 
pect that today's planetesimal belts had too little 
solid mass to overcome the growth slowdown inher- 
ited from their wide orbits. Excepting cases where 
observations conclusively demonstate that a partic- 
ular debris disks and giant planet dynamically inter- 
act (Quillen [2006], Lovis et al. [2006], Kalas et al. 
[2008] ) , the ensemble of debris disks orbiting known 
radial- velocity planet hosts provides an upper limit 
to the semi-major axis of the locations where most 
of the known planets could have formed. 

In f}2]we describe our sample of target stars. In §3] 
and §4j we detail our procedure for reducing Spitzer 
IRS data and detecting infrared excesses. In SjS]we 
discuss our debris disk detections. We focus specif- 
ically on the luminosity and location of the dust in 
f}6] Finally, we list our conclusions in [J71 

2. SAMPLE SELECTION AND OBSERVATIONS 

The main sample of extrasolar planetary systems 
was selected for Spitzer program 40096 (PI: J. Car- 
penter) from the compilation r eported in the Extra- 
solar Planets Encyclopaedia (http://exoplanet.eu) 
as of February 16, 2007. At the time, there were 
182 known extrasolar systems containing 212 plan- 
ets. From this parent sample, we selected extrasolar 
planetary systems that met the following criteria: 

1. The planets were discovered by radial veloc- 
ity techniques, or from transiting surveys and 
later confirmed by radial velocity measure- 
ments. We excluded planets identified from 
microlensing surveys, pulsar timing experi- 
ments, and direct imaging of star forming re- 
gions. 

2. Target stars had distances available from Hip- 
parcos parallaxes. 

3. Targets had either main-sequence or sub- 
dwarf luminosity class to focus the survey of 
solar- analogs. 

6 Based on histograms provided by the Exoplanet Ency- 
clopedia, |http: / / exoplanet . eu| 



4. Targets were bright enough to reach a signal 
to noise ratio of 20 in a spectrophotometric 
bandpass between 30-34 /xm within 20 AOR 
cycles of 120 seconds. 

These criteria narrowed the list to 143 planet 
hosts, all within a distance of 100 pc, containing 
171 planets. Stellar age estimates for the sample 
are available from the detailed spectroscopic work 
by Valenti & Fischer (2005), and/or Ca II H and 
K stellar activity measurements (Saffe et al. [2005], 
Wright et al. [2004]). We cross-correlated this list 
of stars with the Reserve Object Catalog using a 
5x5 arcminute search area and identified 34 extra- 
solar planet hosts with previous IRS observations, 
leaving 109 stars for our IRS observing program. 
We excluded seven stars observed in program 40096 
from our final sample because of poor sky subtrac- 
tion or stray light from nearby sources. 

After the observations were complete, the exis- 
tence of the planet orbiting HD 188753 (Konacki 
2005) was challenged by Eggenberger et al. (2007). 
Here we present our observations of HD 188753 but 
classify it as a non-host. Information about our 
sample, including the number of planets orbiting 
each star and the 70 fim dust luminosity, is given in 
Table 1. Effective temperatures were taken from the 
following sources in order of preference: (1) Valenti 
& Fischer (2005, abbreviated VF05 in Table 1), (2) 
Santos et al. (2004b), and (3) the planet discovery 
paper. 

We also present observations of nine solar-type 
stars with previously unpublished data from other 
IRS programs. From IRS program 41 (PI: G. 
Rieke), we included the planet hosts HD 50554, 
HD 52265, HD 117176 (70 Vir), and HD 134987, 
bringing the total number of planet hosts in our 
sample to 105. Finally, we included the stars with- 
out planets HD 166, HD 33262 and HD 33636 from 
program 41, and HD 105211 and HD 219482 from 
program 2343 (PI: C. Beichman). Our final sample 
includes 111 stars, six of which have no detected 
planets. 

Of the stars in our sample, 109 are FGK dwarfs 
and two, GJ 876 and GJ 581, are M dwarfs. All but 
three of the stars have 70 (im photometry cither 
from the MIPS component of our observing cam- 
paign, described by Bryden et al. (2009), or from 
previous surveys. All but six of our targets are older 
than 1 Gyr, which reduces any spurious statistical 
effects that might result from the well-know correla- 
tion between debris disk brightness and stellar age 
(Su et al. 2006). Thus, our discovered debris disks 
contain information about the planetesimal popula- 
tions of mature planetary systems such as our own. 

The Infrared Spectrograph (IRS) aboard the 
Spitzer Space Telescope observed the 102 targets 
from program 40096 between July 2007 and Jan- 
uary 2008 (IRS campaigns 42-47). For each star, 
we recorded five six-second exposures in the SL1 
module to characterize the stellar photosphere be- 
tween 7 and 14 /im. All stars received a minimum 
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exposure time of 30 seconds (five six-second cycles) 
in the LL2 and LL1 modules, with extra cycles 
and/or ramp time added as needed for faint tar- 
gets to detect the photosphere at the 2-3er level in 
each resolution element. We imposed a maximum 
integration time of 120 sec with 20 cycles to keep 
the AORs to a reasonable length. We also required 
a minimum of five cycles to achieve the redundancy 
needed to identify cosmic rays and bad pixels, and 
to compute uncertainties based on repeatability of 
the spectral extraction. The IRS follow-up obser- 
vations of 70 /im excesses from IRS programs 41 
and 2343 went deeper in the LL1 module, reaching 
S/N ~ 20 per resolution element between 30 and 
34 /im. 

3. DATA ANALYSIS 

Bad pixels were identified using the campaign- 
specific rogue pixel masks provided by the Spitzer 
Science Center. We used the IRSCLEAN package 
to remove bad pixels from the Spitzer pipeline's ba- 
sic calibrated data. Sky subtraction and spectral 
extraction were performed with the SMART data 
reduction package developed by the IRS instrument 
team at Cornell (Higdon et al. 2004). Each star 
was observed at two positions along the slit, Nod 1 
and Nod 2. For each nod, we constructed sky and 
stray-light correction frames by taking the median 
of all exposures performed in the same module but 
opposite nod. 

To make the end-to-end spectrum, we clipped un- 
reliable data points from the beginning and end of 
each order. We retained SL1 data from 7-14/xm, 
LL2 data from 14-20/im, and LL1 data from 20- 
35/im. We imposed a requirement that the final 
spectrum be continuous over order crossings to cor- 
rect for the different slit losses in each order. A 
smooth spectrum was produced by correcting the 
LL2 spectrum with a constant scaling factor to 
match the SL1 flux in the overlap region at 14/im, 
then scaling LL1 fluxes to match LL2 in the 20/im 
overlap region. 

The final step in our data reduction process was 
to create a model of the star's photosphere. For 
each FGK star, we used the Atlas9 grid of model 
atmospheres for stars with solar abundance ratios 
(Kurucz 1992). Following the recommendations of 
Bertone et al. (2004), we used the NextGen model 
atmospheres (Hauschildt et al. 1999) for photo- 
spheres of the M Dwarfs GJ 876 and GJ 3021, which 
have effective temperatures less than 4000 K. 

Photosphere models were fitted using Hipparcos 
parallaxes and Hipparcos/Tycho BT and VT mag- 
nitudes (Perryman et al. 1997; H0g et al. 2000) 
transformed to the Johnson photometric system, 
2MASS J, H and Ks magnitudes (Cutri et al. 
2003), and the measured temperature of each star 
from planet-search spectra. Sources of T e s mea- 
surements are listed in Table 1. When available, 
R and / photometry from Bessell (1990) with as- 
sumed 0.1 mag uncertainties were added to fur- 



ther constrain the photosphere models. Each model 
photosphere was extrapolated to mid-infrared wave- 
lengths using a simple blackbody extension to the 
optical and near- infrared fit. 

In the next section, we describe our method for 
detecting non-photospheric emission, the signature 
of circumstellar dust. 

4. DETECTING EXCESS EMISSION 

To construct the infrared excess spectra, we sub- 
tracted the model photosphere from the IRS spec- 
trum of each star and calculated the fractional flux 
in each resolution element relative to the photo- 
sphere: 

F v (observed) — F v (photosphere) 
^(photosphere) 

Spitzer MIPS 24 /im observations of a sample of 
FGK dwarfs with median age 4 Gyr have demon- 
strated that only ~ 1% of main-sequence Solar- type 
stars have dust with luminosity L& ust / L* > 10 -5 at 
24 /im (Bryden et al. 2006), which originates from 
asteroid-belt analogs at 3-4 AU from Sunlike host 
stars. The fractional excess spectrum of a solar- 
type star with a typical debris disk should stay near 
zero until between 25 and 30 /im, at which point the 
dust gives it a rising slope. With the low likelihood 
of detecting short-wavelength excesses in mind, our 
first step in calculating the fractional excess spec- 
trum was to assume the first 10 resolution elements 
(20.93-22.45 /im) of our clipped LL1 spectra defined 
the photosphere, such that 

i />22.45/jm 

-T- / Udv = 0, (2) 

^ v J20.93fj,m 

where is the frequency interval between 20.93 
and 22.45 /im. 

In the previous section we described using the 
order overlap regions to correct for differential slit 
losses, a procedure that produces a smooth but not 
flux-calibrated spectrum. In order to perform spec- 
trophotometry of our sources it is essential that we 
avoid the unreliable parts of our spectra — the well- 
known SL1 order curvature and the LL1 24 /im flux 
deficit 1 . Pinning the short- wavelength end of LL1 
to the photosphere model defines a flux calibration 
using an instrumentally reliable part of our spectra 
that should contain only photospheric emission. In 
two cases where we did observe a rising fractional 
excess slope shortward of the 20.93 /im LL1 bound- 
ary (HD 166 and HD 219482), we used the first 10 
resolution elements of LL2 (14.24-15.00 /im) to de- 
fine the photosphere. 

Following the approach of Beichman et al. (2006b) 
and Lawler et al. (2009), we defined a filter span- 
ning 30-34 /im which is primarily sensitive to dust 
with T « 115 K. Such dust would be located be- 
tween 8 and 10 AU from a Sunlike star assuming 

1 See a description of these problems at 
http://ssc.spitzer.caltech.edu/irs/features 
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blackbody grains. To search for debris disks, we 
calculate the frequency-weighted mean excess rela- 
tive to the photosphere in the 30-34 /im filter. The 
excess measurement error a based on uncertainties 
propagated through the Spitzer pipeline and photo- 
spheric fitting error is 

i />34/im 

In equation [3l ov,obs is the uncertainty in each reso- 
lution element from the Spitzer pipeline propagated 
through the SMART data reduction, v is the fre- 
quency of each resolution element, and Af is the 
frequency difference between 30 and 34 /im. The 
star flux in equation [3] is simply the frequency- 
weighted average flux of the photosphere model be- 
tween 30 and 34 /im: 

= — — / F v (photosphere) dv. (4) 

Following Beichman et al. (2006b), we impose a 
minimum fractional photosphere fitting uncertainty 
of 2%, which gives the (0.02F*) 2 component of the 
variance. In writing equation [3] we assume that the 
observational errors in each resolution element are 
independent. 

If the errors reported by the Spitzer 
pipeline/SMART and the 2% photosphere-fitting 
uncertainty account for all of the dispersion in our 
infrared excess measurements, we expect only one 
source to have more than a 2<r flux deficit. However, 
a closer look at the distribution of excesses in our 
sample reveals seven stars with 2-3 a flux deficits 
between 30 and 34 /im, indicating that there is 
an extra source of error in our infrared excess 
measurements. Before conclusively identifying any 
debris disks, we follow Beichman et al. (2006b) and 
require that they have a significant infrared excess 
relative to the dispersion of the entire sample, 
not just internal measurement and fitting errors. 
Accordingly, we rank each debris disk candidate 
by the significance of its infrared excess against 
the internal error a. With such a ranking system 
we can identify low-luminosity debris disks in high 
signal-to-noise observations while rejecting debris 
disk candidates that appear to be bright but also 
have large measurement errors. 

We construct a histogram in the variable 

N = Fobs - F * , (5) 
a 

where N is the number of standard deviations a 
from a perfect photosphere in the 30-34 \xm filter 
and -F bs is the frequency- weighted average observed 
flux. We can calculate F ODS by replacing the pre- 
dicted photospheric flux F v (photosphere) in equa- 
tion @] with the observed flux F v (observed) . Figure 
Q] shows the resulting histogram along with its best- 
fit Gaussian function. Since a Gaussian function is 
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Fig. 1. — Debris disks were identified by fitting a Gaus- 
sian distribution, plotted in red, to a histogram of fractional 
30-34 fim excesses. The histogram variable, shown in black, 
is the number of standard deviations <r of the measured ex- 
cess from a perfect photosphere, where a is the frequency- 
weighted quadrature sum of the internal measurement error 
and photosphere fitting error. Our sample has an overall frac- 
tional excess dispersion of £ = 1.93<r (see equation 5). The 
blue cross-hatched region contains our debris disks, which 
have fractional infrared excesses greater than 3E. Our sam- 
ple includes 11 stars with statistically significant 30-34 
excesses. 

an excellent fit to the histogram in Figure [TJ the ex- 
tra error source is random rather than systematic. 
From the standard deviation of our model Gaussian, 
we find that the dispersion of N across our sample 
is £ = 1. 93a. Accordingly, we increase our spec- 
trophotometric error estimates by a factor of 1.93 
to select only debris disks with 3S or greater ex- 
cesses. The cross-hatching in Figure [1] demarcates 
the part of the N distribution where the debris disks 
lie. With the aforementioned detection procedure, 
we find statistically significant 32 /im excesses in the 
spectra of 11 stars. In the next section we discuss 
the spectra of our debris disks. 

5. DEBRIS DISKS 

Table 2 shows the results of our 30-34 /im spec- 
trophotometry. The HD numbers of debris disks 
appear in bold in the table. We detect two new de- 
bris disks orbiting the planet hosts HD 108874 and 
HD 130322. We also measure 30-34 /im emission 
from nine other debris disks previously detected at 
70 /im, four of which orbit planet hosts. The ta- 
ble also gives fractional dust fluxes at 30-34 /im and 
total luminosities for each star in the sample. As- 
suming the dust is a single-temperature blackbody, 
its luminosity is 

-kdust _ F dust kT^ ust (exp [hv/kT dust ] - 1) 
U ~ F* hvTl ' [ ' 

where i*d us t = Fobs — F*. When the dust temper- 
ature is not known, one can use equation [5] to cal- 
culate the minimum dust luminosity, assuming the 
emission peak of the dust is centered in the spec- 
trophotometric filter. For peak dust luminosity at 
32 /im, Xd ust = 115 K. 

One can calculate the true luminosity ratio by 
detecting the debris disk at multiple wavelengths. 
We write equation [6] once each for our 32 /im spec- 
trophotometric fluxes and the MIPS 70 /im values 
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from the literature and solve the resulting system for 
Tdust and Ldust/L*- Note that equation [6] assumes 
that the dust has a single temperature — debris disks 
containing a large range of grain sizes or occupying 
a wide range of distances from the star will have a 
range of dust temperatures. 

Table 2 contains three different types of measure- 
ments of Ldust/i*: true values, lower limits, and 
pseudo-upper limits. For stars without any mea- 
sured infrared excess, we use equation [5] and set 
-Fdust = 3S to calculate the maximum dust luminos- 
ity assuming the excess emission comes from single- 
temperature dust at 115 K. For cold, faint dust un- 
detected by both MIPS and IRS, the true luminosity 
may be higher than the maximum value we calcu- 
lated if the dust emission peaks at very long wave- 
lengths. Lower limits to the dust luminosity arise 
for debris disks with a detection at only one wave- 
length and true values come from detections at mul- 
tiple wavelengths (see S}6]). 

Figures [2] and [3] show the spectra of the debris 
disks in our sample. Each debris disk has little or 
no fractional infrared excess in the IRS SL1 and 
LL2 modules but a positive fractional excess slope 
in the long wavelengths of the LL1 module. Two 
of the debris disks, HD 219482 and HD 166, show 
strong evidence of 24 /jm emission. For compari- 
son, we also show the spectra of three stars without 
any significant infrared excess in Figure These 
spectra tightly hug the line of zero fractional ex- 
cess at all wavelengths and do not show a positive 
fractional excess slope in LL1. Note that the spec- 
trum of HD 50499 shows the SL1 order curvature 
between 10 and 13 /im. As explained in 21 we chose 
the first ten undipped resolution elements of LL1 
(20.93-22.45 /im) to define our photosphere in or- 
der to avoid this order droop. 

Our overall detection rate of debris disks is 11%, 
but the non-uniform sensitivity and the inclusion 
of 70 /im follow-up observations in our survey bias 
this statistic. Because of the non-uniformity of our 
observations, we do not attempt to calculate an er- 
ror bar for our detection rate. However, our results 
agree with the IRS study of Beichman et al. (2006b) , 
who found 32 /im excesses around 12 ± 5% of stars 
in a sample that also included follow-ups of 70 /im 
excesses. Similarly, Lawler et al. (2009) report a 
32 /im excess detection rate of 11.8 ± 2.4%. 

The most intrinsically luminous debris disk in our 
sample is that orbiting the planet host HD 202206, 
which has L^ust/L* > 10~ 4 . For comparison, the 
well-studied debris disks orbiting the planet hosts 
HD 69830 and HD 82943 have luminosity ratios 
Ldust/L* of 2 x 10~ 4 and 1.2 x 10 -4 respectively 
(Beichman et al. 2005). All but two of the 111 stars 
in our sample have a minimum 3E detectability of 
Ld-ast/L* < 10 -4 ; thus only ~ 1% of the sample has 
a fractional dust luminosity brighter than lO -4 !/*. 
Our result agrees with the 2 ± 2% detection rate 
for excesses of I/dust/^* calculated by Bryden et al. 
(2006). 
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Fig. 2. — The debris disks in our sample have fractional ex- 
cess spectra that rise toward longer IRS wavelengths. Only 
a few have any noticeable rise near 24 /im, consistent with 
the paucity of 24 /an excesses observed by the MIPS in- 
strument. Several spectra show the well-known SL1 order 
droop. HD 108874, HD 117176, HD 130322 and HD 202206 
are planet hosts. 
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Fig. 3. — Fractional excess spectra of detected debris disks, 
continued. No SL data points are included for HD 33262 
because the SL1 spectra were partially saturated. HD 50554 
and HD 52265 are planet hosts. 



6. DUST TEMPERATURE AND LOCATION 

We will now determine the location of the debris 
for the sample of 11 stars with infrared excesses de- 
tected at 32 /mi. Of our 11 detected debris disks, 
nine also have excess emission at 70 /im (Beichman 
et al. [2006a], Trilling et al. [2008], Bryden et al. 
[2009], Kospal et al. [2009]). Table 2 lists the black- 
body temperature of each debris disk detected at 
multiple wavelengths. Beichman et al. (2006a) de- 
tected HD 219482's dust at both 24 and 70 (im and 
calculated Td us t — 82 ± 3 K for large grains, in ex- 
cellent agreement with our value of 81 ± 3 K. Our 
79 ± 3 K temperature for the dust orbiting HD 166 



G 



Dodson-Robinson et al. 



HD50499 (no excess) 



HD142 [no excess) 



15 20 25 30 35 40 

Wavelength (micron) 



Fig. 4. — Fractional excess spectra of stars with no signifi- 
cant 30-34 fim emission. These stars show no noticeable rise 
in excess toward long wavelengths. 

does not agree with Tanner et al. (2009), who cal- 
culated a lower limit of 104 K using the previously 
published 70 pm detection (Trilling et al. 2008) and 
the 160 pm upper limit. However, grain tempera- 
ture varies with size roughly as 
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where r gr is the grain radius and a is the semima- 
jor axis of the debris orbit (Krugel 2003). There- 
fore each photometric band is most sensitive to a 
certain combination of grain size and distance from 
the star. Different pairs of photometric bands may, 
therefore, lead to different grain temperature dis- 
tance estimates (Carpenter et al. 2009). 

By combining our IRS observations with previ- 
ous MIPS data, we have obtained temperature es- 
timates for 15 planet hosts and five non-host stars. 
Figure [5] shows luminosity ratio vs. dust temper- 
ature for all debris disks in our sample, whether 
detected in this study or previous 70 /im MIPS 
work. The warmest temperature of a debris disk de- 
tected at multiple wavelengths is 91 K for the dust 
around HD 33262. Only four stars may possibly 
host dust warmer than 100 K— HD 142, HD 33262, 
HD 108874 and HD 130322 have temperature upper 
limits of 117 K, 103 K, 120 K and 120 K respectively. 
Our finding that debris disks around solar-type stars 
are primarily cold agrees with Hillenbrand et al. 
(2008), who find dust temperatures < 45 — 85 K 
for the 25 debris disks in the FEPS Spitzer Legacy 
Science Program. 70 /jm debris systems with tem- 
perature upper limits only may be true Kuiper Belt 
analogs, as most have Td us t < 70 K. 

With the exception of three of the upper limits 
quoted above, all dust temperatures are below the 
water ice sublimation temperature of T ~ 110 K 
in a zero-pressure medium. However, in optically 
thin disks photodesorption is the dominant volatile 
removal mechanism and is efficient well beyond the 
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Fig. 5. — Dust luminosity vs. temperature for all debris 
disks in our sample. Black circles indicate the debris disks 
detected at both 32 and 70 fim. Red error bars indicate the 
possible temperature range for debris disks detected only at 
32 fim while right-pointing red arrows denote that the lu- 
minosity ratios for these debris disks are lower limits only. 
Purple points with double arrows denote the debris disks 
detected only at 70 /im, which gives upper limits for tem- 
peratures and lower limits for luminosities. Only four stars, 
HD 142, HD 33262, HD 108874 and HD 130322, may have 
debris rings with dust hotter than 100 K (see £|6t . 

thermal ice line (Grigorieva et al. 2007). Indeed, the 
survival timescale for ice deposited in monolayers on 
the surfaces of small grains in an optically thin disk 
is of order a few minutes (Chen et al. 2005). Assum- 
ing, accordingly, that (a) the dust in each system is 
rocky or carbonaceous with an albedo near zero and 
(b) most grains are larger than 10 /jm and so have 
absorption coefficients Q ~ 1 for visible and near- 
IR irradiation, we can find the distance of the dust 
from its host star in the blackbody approximation, 
appropriate for large grains: 



z J dust 



(8) 



where i?* is the star radius. For our solar-type stars 
with 60 < T dust < 100 K, we find 7AU < a < 20AU. 
Large grains in the debris disks we have discovered 
at 32 /im occupy what would be the region between 
Saturn and Uranus in our own Solar System. 

Smaller grains in the 60-100 K temperature range 
would reside on larger orbits. We estimate the min- 
imum grain size in our debris disks from the radia- 
tion pressure blowout limit (Artymowicz 1988) of 



WnGM^cp 1 



(9) 



where r m - m is the minimum grain size, L* and M* 
are the star luminosity and mass, Q pr is the ra- 
diation pressure coupling coefficient, and p is the 
internal density of the grains. Q pr = 1 is a good ap- 
proximation for Solar- type stars (Krugel 2003). For 
silicate grains with p — 3.3 g cm~ 3 , a mm — 0.1 pm 
Icy grains with p = 1.0 g cm 
type star have minimum size a 
equation which applies to grain sizes from ap- 
proximately 0.03 to 30 ^m (Kruegel 2003), we find 
that 70 K silicate grains at the blowout limit orbit 
~ 240 AU from Solar- type stars. Low-density icy 
grains, which would be destroyed quickly by pho- 



3 orbiting a Solar- 
0.3 pm. Using 
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todesorption, are limited to within <~ 142 AU from 
their host stars. 

In the introduction, we emphasized the lack of any 
dynamical connection between short-period planets 
and cold dust — planets detected by radial velocity 
neither migrated through and destroyed our plan- 
etesimal belts, nor dynamically stir them today. 
Based on the 240 AU maximum orbital radius of 
grains with temperatures near 70 K, a representa- 
tive temperature for our detected debris disks, we 
can restrict the formation locations of the known 
radial-velocity planets to within ~ 240 AU. In 
any individual debris system where the bulk of the 
infrared excess emission comes from large, near- 
blackbody grains, we can restrict the giant planet 
formation zone even further, excluding regions be- 
yond 15 AU. In such systems, giant planet forma- 
tion by core accretion would have been restricted 
to a narrow belt near the ice line, consistent with 
the theoretical predictions of Dodson-Robinson et 
al. (2009b). So far, the orbital radii of the 15 debris 
disks orbiting planet hosts in this sample provide 
the only purely observational constraint on where 
the radial-velocity planets could have formed. 

7. CONCLUSIONS 

We searched a sample of 111 stars, including 105 
planet hosts, for 30-34 /zm excess emission indicat- 
ing circumstellar dust. Out of 111 stars, 11 showed 
detectable 32 fim spectrophotometric excess. Our 
debris disk detection rate of 11% is consistent with 
previous IRS and MIPS 70 /xm studies of solar-type 
stars. We have identified two new debris disks orbit- 
ing HD 108874 and HD 130322 and calculated dust 
temperatures and luminosities for the debris disks 
that had been previously detected at 70 itm. 

Our IRS data show that the majority of debris 
systems contain cold dust with Td us t < 100 K. 
Here we confirm earlier results from 24 /im pho- 
tometry (Bryden et al. 2006) and IRS spectroscopy 
(Hillenbrand et al. 2008) demonstrating that ex- 
cess emission is usually detectable only at wave- 
lengths longer than 30 fim. The exception in our 
survey that has noticeable excess at 20 xtm, HD 166, 
is cooler than 80 K. Lower limits to the circum- 
stellar dust temperatures of Solar-type stars come 
from the MIPS 160 /xm photometry of Tanner et 
al. (2009); they find that T dust > 27 K for the ma- 
jority of the debris disks in their sample. While 
all of our debris disks are colder than the ice sub- 
limation temperature, we expect them to contain 
mainly bare rock/carbonaceous grains since photo- 
evaporation of ice is extremely efficient in optically 
thin disks (Grigorieva et al. 2007). 

The cold temperatures of most circumstellar dust 
disks make them excellent candidates for follow- 
up with the Herschel Space Observatory (Pilbratt 
et al. 2010). The PACS medium-resolution spec- 
trometer, which covers 55-210 /im at a resolution of 
R = A/AA =1000-4000 and also includes photomet- 
ric bands, will cover the peak blackbody emission of 



such disks and provide detailed dust spectral energy 
distributions. Circumstellar dust SEDs constrain 
the location of the emitting dust and show whether 
it occupies a narrow ring — the simplifying assump- 
tion used in this work — or a broad, multi-AU band 
(e.g. Hillenbrand et al. 2008). If any debris disks 
contain icy dust produced by very recent collisions, 
PACS may be able to detect H2O gas produced by 
photoe vapor ation of the grain mantles. 

One question open for debate is why the planetes- 
imal belts surrounding our host stars did not accu- 
mulate into giant planets. While planets may shape 
debris disk edges (Quillen 2006), they will eject 
planetesimals within a few Roche lobe radii of their 
orbits (Bryden et al. 2000), so planets and circum- 
stellar dust should not coincide. Planet-formation 
simulations show that formation of ice giants such 
as Uranus and Neptune is robust out to ~ 15 AU 
in disks of several times the minimum-mass solar 
nebula (Dodson-Robinson et al. [2009a], Dodson- 
Robinson & Bodcnhcimer [2010]). Blackbody grains 
with Tdust = 70 K reside only 15 AU from solar- type 
stars. 

If most of our debris disk luminosity comes from 
blackbody grains, a simple explanation for the pres- 
ence of planetesimal belts in locations that could 
theoretically host planets is that the majority of 
protostellar disks do not have enough mass to form 
ice giants near ~ 15 AU. Gas and ice giant forma- 
tion by core accretion would then be restricted to 
a narrow region around the ice line in most disks, 
as predicted by Dodson-Robinson et al. (2009b). 
The Sun's ability to form Uranus and Neptune 
in the trans-Saturnian region (Dodson-Robinson & 
Bodenheimer 2010) would be the unusual result of 
an abnormally massive solar nebula. 

If most of the infrared excess emission is produced 
by small grains near the radiation-pressure blowout 
limit, we can in general confine the ensemble of 
radial-velocity planets to formation zones less than 
240 AU from host stars. We have therefore pre- 
sented the first fully observational constraints on 
the formation locations of the ensemble of radial- 
velocity planets. 
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TABLE 1 

Observed Stellar Sample 



HD 


Name 


Spectral Type 


T cS (K) 


T e g Reference 


V (mag) 


Dist. (pc) 


Planets 


%ist x io 5 (70 /urn) 
* 


70 Reference 




GJ 581 


M3 


3310 


Bessell (1995) 


10.57 


6.3 


3 


8.49 


Kospal ct al. (2009) 




GJ 876 


M4 


3130 


Bessell (1995) 


10.17 


4.7 


3 


<3.2 


Kospal ct al. (2009) 




Hip 14810 


G5 


5485 


Wright ct al. (2009) 


8.52 


52.9 


3 




142 




F7V 


6248 


VF05 a 


5.70 


25.6 


1 


0.8 


B09 b 


166 




KOV 


5577 


VF05 


6.13 


13.7 





6.0 


Trilling et al. (2008) 


1237 


GJ 3021 


G8.5Vk: 


5360 


Santos et al. (2004b) 


6.59 


17.6 


1 


<1.0 


B09 


3651 




KOV 


5220 


VF05 


5.80 


11.1 


1 


<1.1 


B09 


4203 




G5 


5701 


VF05 


8.69 


77.8 


1 


<23.9 


B09 


4208 




G7VFe-lH-05 


5600 


VF05 


7.79 


32.7 


1 


<2.9 


B09 


8574 




F8 


6049 


VF05 


7.11 


44.2 


1 


<1.5 


B09 
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109 Psc 


G5IV 
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K1V 


5150 


VF05 


6.17 
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1 
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G5IV 
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B09 
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Rocha-Pinto & Maciel (1998) 
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B09 
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6188 


VF05 


5.06 
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GO 
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Santos et al. (2004b) 
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Trilling et al. (2008) 
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Trilling et al. (2008) 
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Kospal ct al. (2009) 
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HD 41004 A 
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Santos et al. (2002) 


8.65 


43.0 


1 


<6.9 


B09 


45350 




G5 


5616 


VF05 


7.88 


48.9 


1 


<7.2 


B09 


46375 




K1IV 


5285 


VF05 


7.84 


33.4 


1 


17.8 


Kospal et al. (2009) 


49674 




GO 


5662 


VF05 


8.10 


40.7 


1 


<4.5 


B09 


50499 




G1V 


6069 


VF05 


7.21 


47.3 


1 


1.6 


Kospal ct al. (2009) 


50554 




F8V 


5928 


VF05 


6.84 


31.0 


1 


5.85 


Trilling et al. (2008) 


52265 




GOV 


6076 


VF05 


6.30 


28.1 


1 


2.69 


Trilling et al. (2008) 


63454 




K3V 


4841 


VF05 


9.37 


35.8 


1 


<16.0 


B09 


65216 




G5V 


5666 


Santos et al. (2003) 


7.97 


35.6 


1 


<4.4 


B09 


68988 




GO 


5960 


VF05 


8.20 


58.8 


1 


<4.2 


B09 


70642 




G6VCN+05 


5705 


VF05 


7.18 


28.8 


1 


<6.3 


B09 


72659 




GO 


5919 


VF05 


7.46 


51.4 


1 


<2.2 


B09 


73256 




G8IV-VFe+05 


5570 


Udry et al. (2003) 


8.08 


36.5 


1 


<5.5 


B09 


74156 




GO 


6067 


VF05 


7.61 


64.6 


3 


<3.4 


B09 


75289 




F9VFe+03 


6095 


VF05 


6.36 


28.9 


1 


<14.9 


B09 


76700 




G6V 


5668 


VF05 


8.13 


59.7 


1 


<3.7 


B09 


81040 




GO 


5700 


Sozzctti ct al. (2006) 


7.72 


32.6 


1 


<3.2 


B09 



TABLE 1 

Observed Stellar Sample 



83443 




KOV 


5453 


VF05 


8.24 


43.5 


1 


<9.0 


B09 


88133 




G5IV 


5494 


Fischer ct al. (2005) 


8.06 


74.5 


1 






89307 




GOV 


5897 


VF05 


7.01 


30.9 


1 


<2.8 


B09 


89744 




F7V 


6291 


VF05 


5.74 


39.0 


1 


<0.9 


B09 


93083 




K2IV-V 


4995 


Levis et al (2005") 


8.33 


28.9 


1 


<9.2 


B09 


99492 




K2V 


4954 


VF05 


7.57 


18.0 


1 


<5.1 


B09 


102117 




G6V 


5695 


VF05 


7.47 


42.0 


1 


<16.8 


B09 


102195 




KO 


5291 


Melo et al. (2007) 


8.06 


29.0 


1 


<9.2 


B09 


105211 




F2V 


6600 c 


Peletier (4989*1 


4.15 


19.7 





6.9 


Beichman et al (2006a "l 


107148 




G5 


5797 


VF05 


8.01 


51.3 


1 


<5.6 


B09 


108147 




F8VH+04 


6156 


VF05 


7.00 


38.6 


1 


<23.9 


B09 


108874 




G5 


5550 


VF05 


8.74 


68.5 


2 


<13.6 


B09 


109749 




G3V 


5903 


Fischer et al. (2006) 


8.08 


59.0 


1 


<6.2 


B09 


111232 




G8VFe-10 


5494 


Santos ct al. (2004b) 


7.61 


28.9 


1 


<6.1 


B09 


114386 




K3V 


4819 


VF05 


8.73 


28.0 


1 


<6.8 


B09 


114762 




F9V 


5952 


VF05 


7.30 


40.6 


1 


<3.3 


B09 


114783 




KO 


5135 


VF05 


7.56 


20.4 


1 


<3.0 


B09 


117176 


70 Vir 


G5V 


5544 


VF05 


5.00 


18.1 


1 


0.97 


Trilling et al. (2008) 


117207 




G7IV-V 


5723 


VF05 


7.26 


33.0 


1 


<2.1 


B09 


117618 




GOV 


5963 


VF05 


7.17 


38.0 


1 


<2.3 


B09 


118203 




KO 


5695 


da Silva et al. (2006) 


8.05 


88.6 


1 


<13.2 


B09 


130322 




KOIII 


5308 


VF05 


8.04 


29.8 


1 


<7.9 


B09 


134987 


23 Lib 


G5V 


5750 


VF05 


6.45 


25.7 


2 


<2.70 


Trilling et al. (2008) 


136118 




F8 


6097 


VF05 


6.93 


52.3 


1 


<1.3 


B09 


142415 




G1V 


5901 


VF05 


7.34 


34.6 


1 


<10.9 


B09 


143761 


p Cr B 


GOV 


5823 


VF05 


5.40 


17.4 


1 


<0.6 


B09 


145675 


14 Her 


KOV 


5347 


VF05 


6.67 


18.1 


2 


<1.0 


B09 


147513 




G1VH-04 


5929 


VF05 


5.38 


12.9 


1 


<1.6 


B09 


149026 




GOIV 


6147 


Sato et al. (2005) 


8.15 


78.9 


1 


<5.2 


B09 


149143 




GO 


5884 


Fischer et al. (2006) 


7.89 


63.5 


1 


<4.6 


B09 


154857 




G5V 


5605 


VF05 


7.24 


68.5 


1 


<4.7 


B09 


159868 




G5V 


5623 


VF05 


7.24 


52.7 


1 


<9.9 


B09 


160691 


p Arae 


G3IV-V 


5784 


Santos et al. (2004a) 


5.15 


15.3 


4 


<0.9 


B09 


162020 




K2V 


4844 


VF05 


9.10 


31.3 


1 


<28.6 


B09 


168443 




G6V 


5579 


VF05 


6.92 


37.9 


2 


<14.7 


B09 


169830 




F7V 


6221 


VF05 


5.91 


36.3 


2 


<1.5 


B09 


177830 




KO 


4948 


VF05 


7.18 


59.0 


1 


<1.4 


B09 


183263 




G2IV 


5936 


VF05 


7.86 


52.8 


2 


<9.1 


B09 


185269 




GOIV 


5980 


Johnson et al. (2006) 


6.68 


47.4 


1 


<2.4 


B09 


186427 


16 Cyg B 


G3V 


5674 


VF05 


6.20 


21.4 


1 


<1.3 


B09 


187085 


GOV 


6075 


VF05 


7.22 


45.0 


1 


2.3 


Kospal ct al. (2009) 


187123 




G5 


5814 


VF05 


7.83 


47.9 


2 


<10.6 


B09 


188015 




G5IV 


5745 


VF05 


8.24 


52.6 


1 


<29.4 


B09 


188753 


HD 188753 A 


G8V 


5750 


Konacki (2005) 


7.40 


44.8 


d 


<8.7 


B09 


190360 




G7IV-V 


5551 


VF05 


5.71 


15.9 


2 


<9.7 


B09 


192263 




K2V 


4975 


VF05 


7.79 


19.9 


1 


5.4 


B09 


192699 




G5 


5220 


Johnson ct al. (2007) 


6.45 


67.4 


1 


<1.1 


B09 


195019 




G3IV-V 


5788 


VF05 


6.91 


37.4 


1 


<2.1 


B09 


196050 




G3V 


5892 


VF05 


7.50 


46.9 


1 


<2.7 


B09 


196885 




F8IV 


6185 


VF05 


6.40 


33.0 


1 


<1.0 


B09 


202206 




G6V 


5787 


VF05 


8.08 


46.3 


2 


14.3 


B09 


208487 




G2V 


6067 


VF05 


7.47 


44.0 


1 


<4.9 


B09 



TABLE 1 

Observed Stellar Sample 



210277 


GOV 


5555 


VF05 


212301 


F8V 


6256 


Lo Curto ct al. (2006) 


213240 


G0/G1V 


5967 


VF05 


216435 


GOV 


5999 


VF05 


216770 


G9VCN+1 


5229 


Mayor et al. (2004) 


217107 


G8IV 


5704 


VF05 


219482 


F6V 


6240 


Flower (1996) 


222404 7 Cep 


K1IV 


4916 


Santos ct al. (2004b) 


222582 


G5 


5726 


VF05 


224693 


G2V 


6037 


Butler et al. (2006) 



VF05 = Valcnti & Fischer (2005) 
BOO = Brydcn ct al. (2009) 

Based on tabulated values for vs. V — K 

Planet's existence under debate; sec Eggcnbcrgcr ct al. (2007b) 



6.63 


21.3 


1 


<0.8 


B09 


7.76 


52.7 


1 


<4.9 


B09 


6.80 


40.7 


1 


<1.1 


B09 


6.03 


33.3 


1 


2.2 


B09 


8.10 


37.9 


1 


<19.0 


B09 


6.18 


19.7 


2 


<1.5 


B09 


5.66 


20.6 





2.8 


Bcichman et al. 


3.23 


13.8 


1 


<0.3 


B09 


7.68 


41.9 


1 


<2.3 


B09 


8.23 


94.1 


1 


<8.0 


B09 



(2006a) 



Spectra of Debris Disks 
TABLE 2 

IRS 30-34 /j,m Spectrophotometric Results 



HD 


Name 


F, 


-^obs 


£ 


-Fdust/ i 7 * 


_dusjL X 10° a 


T dust (K) b 






16.1 


17.1 


0.99 


0.06 


> 12 


< 83 




GJ 876 


24.8 


25.3 


0.67 


0.02 


< 5.5 






Hin 1/L&10 


7.5 


7.0 


0.93 


-0.07 


\ O. I 




142 




30.9 


30.8 


2.67 


0.00 


> 3.0 


< 119 


166 




78.9 


119 


2.00 


0.50 


7.4 


79 ± 3 


1237 


OT 3021 


39.9 


38.7 


1.71 


-0.03 


< 1.7 




3651 




98.6 


101.0 


2.29 


0.02 


< 0.8 




4203 




5.8 


5.9 


0.40 


0.01 


< 2.0 




4208 




12.8 


13.1 


0.87 


0.02 


< 2.0 




8574 




19.2 


17.7 


1.57 


-0.08 


<r 1 Q 




1 0RQ7 


1 HQ P«r> 
iuy r 




48 1 


1 Q8 


01 
-U.U1 


/ 1 o 




11964 




8.3 


8.3 


1.39 


0.00 






12661 




18.7 


19.6 


0.86 


0.05 


<r 1 A 




13445 


GJ 86 


90.3 


88.3 


3.00 


-0.02 


< 1.3 




16141 




32.4 


31.4 


1.72 


-0.03 


<r 1 A 




17051 


J-l-LV. OIU 


107.4 


108.9 


2.34 


0.01 


< 0.5 




19994 




126.0 


124.0 


2.91 


-0.02 


> 6.9 


< 83 


20367 




76.0 


76.1 


1.43 


0.00 


< 0.4 




20782 




20.3 


20.0 


1.29 


-0.02 


<r 1 7 




23079 




20.1 


21.0 


1.11 


0.05 


\ -L.O 




23127 




5.5 


4.9 


0.65 


-0.10 


< 2.8 








24.9 


25.0 


1.16 


0.00 


<r 1 1 




27442 


e Ret 


241.0 


240.9 


3.52 


-0.30 


< 0.4 




27894 




2.5 


2.5 


0.67 


-0.02 


< 12 




28185 




13.2 


13.1 


0.59 


-0.01 


\ -L.O 




30177 




8.4 


8.0 


0.51 


-0.05 


\ -L.O 






C Dor 


177.3 


192.9 


1.64 


0.09 


0.8 


91 ± 12 


33283 




8.5 


9.5 


0.81 


0.12 


< 2.2 




33564 




109.6 


112.3 


2.78 


0.03 


< 0.5 








24.4 


29.7 


0.85 


0.22 


5.2 


61 i 3 


37124 




16.4 


14.9 


1.18 


-0.09 


< 2.2 




37605 




7.6 


8.4 


0.51 


0.11 


< 2.3 




40979 




27.4 


28.7 


0.90 


0.05 


^> -L.U 


67 


41004 


HD 41004 A 


11.7 


12.2 


0.61 


0.04 


< 2.5 




*±OOOVJ 




u.^ 




47 


01 
-U.U1 


< 2 2 




1UO ( o 




10 7 


_L _L .O 


u.uu 


07 
U.U i 




s 74 


49674 




10.3 


10.6 


0.40 


0.03 


\ -L 




50499 




18.1 


18.3 


0.65 


0.01 


> 1.8 


< 68 


50554 




25.6 


31.0 


0.68 


0.21 


6.2 


58 ± 3 


522R5 

* > Jmd \J KJ 




41.0 


44.5 


0.66 


0.09 


2.9 


56 ± 4 


63454 




7.1 


7.3 


0.47 


0.03 


< 3.5 




65216 




11.4 


10.1 


0.63 


-0.11 


<r 1 

\ -L.U 




uoyoo 




8.3 


8.9 


0.32 


0.08 


1 
\ ±.u 




70642 




23.5 


23.0 


0.63 


-0.02 


< 0.8 




72659 




16.5 


17.7 


0.67 


0.08 


< 1.0 




73256 




12.3 


11.8 


0.52 


-0.05 






74156 




12.2 


11.5 


0.58 


-0.06 


< 1.1 




75289 




38.1 


31.7 


1.54 


-0.17 


< 0.9 




76700 




10.4 


10.0 


0.70 


-0.05 


< 2.0 




81040 




16.0 


16.5 


0.84 


0.03 


\ -L.U 








10 8 


y.o 


41 


-0 12 


\ -L.O 




OOlOO 




14 7 


14 


5fi 


0^ 

-u.uo 


■\ -L.O 




89307 




25.5 


25.0 


0.84 


-0.02 


< 0.9 




89744 




53.7 


53.8 


2.47 


0.00 


< 1.0 








10.0 


10.0 


0.46 


0.00 


< 2.1 




99492 




30.6 


31.4 


1.08 


0.03 


< 0.8 




102117 




19.5 


20.4 


1.09 


0.05 


<r 1 7 

\ -L . ( 




1 021 95 




13.3 


13.9 


0.42 


0.04 


< 1.1 




105211 


?7 ru 


178.5 


215.6 


3.30 


0.21 


6.7 


51 i 2 


1 071 AR 




10 4 


y.u 


fi8 
u.uo 


OR 

-u.uo 


^ i 7 
\ -L. 1 




108147 




29.4 


29.2 


0.99 


-0.01 


< 0.8 




108874 




5.9 


7.9 


0.33 


0.32 


> 1.2 


57 < T < 120 


109749 




8.0 


7.5 


0.61 


-0.06 


< 1.8 




111232 




54.4 


54.3 


1.26 


0.00 


< 0.7 




114386 




12.5 


13.7 


0.71 


0.09 


< 3.0 




114762 




19.3 


19.7 


1.13 


0.02 


< 1.5 




114783 




21.1 


20.7 


1.39 


-0.02 


< 2.6 




117176 


70 Vir 


173.0 


182.1 


2.87 


0.05 


1.0 


72 ± 10 


117207 




21.1 


18.6 


1.61 


-0.12 


< 2.3 




117618 




21.2 


20.4 


1.04 


-0.04 


< 1.3 




130322 




13.0 


16.1 


0.58 


0.24 


> 7.9 


61 < T < 120 


134987 


23 Lib 


39.9 


40.4 


0.76 


0.01 


< 0.5 
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TABLE 2 

IRS 30-34 fj,M Spectrophotometric Results 



136118 




26.7 


26.3 


0.88 


-0.01 


< 0.8 


142415 




17.6 


19.2 


0.59 


0.09 


< 0.9 


143761 


p Cr B 


87.2 


87.6 


2.19 


0.01 


< 0.6 


145675 


14 Her 


53.2 


50.7 


2.57 


-0.05 


< 1.9 


147513 




110.5 


115.4 


2.78 


0.05 


< 0.7 


149026 




7.8 


7.9 


0.35 


0.02 


< 1.0 


149143 




10.8 


11.5 


0.31 


0.06 


< 0.8 


154857 




29.1 


29.5 


0.94 


0.01 


< 1.0 


159868 




23.71 


21.93 


0.78 


-0.075 


< 0.6 


160691 


fA Arac 


247.1 


247.3 


2.36 


0.00 


< 0.3 


162020 




10.2 


9.5 


0.36 


-0.07 


< 1.9 


168443 




30.4 


29.8 


1.70 


-0.02 


< 1.7 


169830 




50.2 


51.3 


1.31 


0.02 


< 0.5 


177830 




43.4 


42.4 


2.48 


-0.02 


< 2.6 


183263 




9.7 


10.8 


0.51 


0.12 


< 1.2 


185269 




34.5 


35.0 


1.43 


0.02 


< 1.0 


186427 


16 Cyg B 


52.6 


54.2 


3.74 


0.03 


< 2.1 


187085 




16.7 


16.9 


0.80 


0.01 


> 2.3 


187123 




12.5 


12.8 


0.57 


0.03 


< 1.2 


188015 




9.1 


7.5 


0.91 


-0.18 


< 3.0 


188753 


HD 188753 A 


24.6 


24.3 


1.17 


-0.01 


< 1.6 


190360 




91.1 


94.4 


1.83 


0.05 


< 0.6 


192263 




22.4 


21.9 


1.04 


-0.02 


> 5.1 


1 AO/ 1 / 1 / 1 

192699 




84.0 


85.6 


2.48 


0.02 


< 1.2 


195019 




31.2 


32.7 


0.80 


0.03 


< 0.8 


1 ^ 1 /" / 1 "~ /I 

196050 




14.9 


16.4 


1.52 


0.10 


< 2.7 


196885 




25.9 


1.00 


0.00 


< 0.9 




202206 




10.4 


14.3 


0.35 


0.38 


12 


208487 




13.5 


12.3 


0.69 


-0.09 


< 1.2 


210277 




53.9 


55.3 


1.98 


0.03 


< 1.3 


212301 




11.2 


11.7 


0.64 


0.04 


< 1.3 


213240 




34.2 


34.4 


1.50 


0.01 


< 1.2 


216435 




59.3 


67.0 


2.81 


0.13 


> 2.1 


216770 




10.2 


10.4 


0.55 


0.02 


< 1.9 


217107 




56.9 


55.1 


3.64 


-0.03 


< 1.9 


219482 




66.5 


91.0 


1.15 


0.37 


3.7 


222404 


7 Cep 


1622.0 


1637.0 


8.36 


0.01 


< 0.2 


222582 




12.4 


12.1 


0.59 


-0.03 


< 1.3 


224693 




6.3 


6.1 


0.62 


-0.03 


< 2.3 



< 67 

< 63 
57 ±3 

< 70 
81±3 C 



a Upper limits to the dust luminosity arc calculated assuming a single-temperature blackbody with Td us t — 115 K. For cold dust 
undetected by both MIPS and IRS, the true luminosity ratio may be higher than the value quoted here. Sec ^5]for further discussion 
of the relationship between observed 30-34 /j,m fluxes and grain luminosity. 

' Grain temperature estimates are based on the blackbody approximation. For a discussion of temperature as a function of grain 
size, sec §E\ 

C Our temperature measurement for the circumstellar dust of HD 219482 agrees with that of Beichman ct al. (2006), who calculated 
Td us t — 82 ± 3 K based on 24 and 70 fim photometry. 



